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A strong coupling expansion for a two-band Hubbard model on two sites with nearly degenerate states 
is considered. A comparative analysis is performed for different schemes of perturbation theory which 
are applicable to systems with nearly degenerate states. A fourth order approach which builds on a four- 
dimensional low-energy subspace with nearly degenerate states captures accurately the transition from an 
antiferromagnetic to a ferromagnetic ground state at large on-site Coulomb interaction. 
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1 Introduction 

The emergence and antagonism of ferromagnetic and antiferromagnetic ground states in strongly corre- 
lated electron systems has always been considered to be of prime importance for a thorough understanding 
of the magnetic states in narrow band systems such as the transition metal oxides. In a large number of 
cases the nature of the magnetic coupling is successfully predicted by appropriate spin-orbital models. For- 
mally the magnetic coupling results from second order expansion in some hopping amplitude parameter t, 
and therefore applies to systems in the strong coupling regime with strong on-site Coulomb interaction U 
(with respect to the kinetic energy t, i.e., t <C U). In that case spin-orbital models provide a framework to 
analyze a wealth of properties (for a recent reference see (Dies et al. Uj and references therein). Neverthe- 
less, many systems are rather tracked down to be in the intermediate coupling regime, where, for example, 
phase transitions between ferromagnetic and antiferromagnetic phases can take place (see, e. g., (2)). Such 
situations are not covered by the standard spin-orbital models, and higher orders in the expansion in t are 
needed. Unfortunately this often implies to take nearly degenerate states into account. In fact, the determi- 
nation of the magnetic coupling of nearly degenerate states beyond second order perturbation theory has 
been rarely addressed in solid state physics. Our paper aims at filling this gap. 

Whereas small clusters with few electronic states can be exactly diagonalized to characterize the elec- 
tronic state, larger clusters with millions of states require different approaches. Already a cluster of three 
transition metal ions (with six d electrons per site) and two oxygen ions provides more than 600 thousand 
low energy states which are nearly degenerate. To be explicit, this challenge is encountered in the quasi- ID 
Mott-insulator Ca3Co2C>6, when effective interaction parameters are calculated from a minimal cluster of 
three Co sites (only counting the 3d states) and two O ligand sites (exclusively considering 2p states), and 
a microscopic analysis of the dominant ferromagnetic coupling between high-spin next-nearest neighbor 
Co 3+ sites has not yet included all low-energy states. In such a system one would have to implement a per- 
turbation theory of at least 5 th order (due to the necessity to include ring exchange). Moreover, an adequate 
perturbation theory should be fully qualified to cope with a large number of nearly degenerate states in the 
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low energy sector. Such a situation with several or many nearly degenerate low energy states, which are 
separated by a large energy scale U from the high energy sector, is typical for multi-band transition metal 
oxides and the question arises if well established perturbation theories are available with which electronic 
systems with almost or nearly degenerate ground states can be reliably handled Q. 

Schemes for the implementation of a perturbation theory with degenerate states in the low energy sector 
have been known for a long time and are well presented in textbooks on quantum mechanics and review 
articles (e.g., see the Refs. B HTT1 ). The extension to almost degenerate states seems to be straightforward 
with either of two alternative procedures: (i) one formally assigns a degenerate "unperturbed" Hamiltonian 
to those eigenstates which are nearly degenerate 05] [7] . The perturbation expansion is then set up with 
respect to this artificially constructed, unperturbed degenerate Hamiltonian with the perturbation being 
comprised of the terms which lift the degeneracy and the original perturbation terms which connect to the 
high energy sector. This scheme, which introduces an artificial degeneracy, is in general not physically 
motivated and may produce inferior convergence as we will demonstrate in Sec. [3] (ii) Alternatively, 
one first diagonalizes the unperturbed Hamiltonian jointly with that part of the perturbation which can be 
projected onto the low energy sector of the Hilbert space (in which the degenerate and almost degenerate 
states live). The part of the perturbation which connects to the high energy sector can then be treated 
in standard perturbation theory J6]. However this latter scheme is valid only up to third order where 
intermediate states in this scheme are exclusively states from the high energy sector. In fourth order, 
this scheme breaks down because intermediate states are then in turn allowed to be taken from the low 
energy sector. We will shortly introduce this procedure in the following section. The standard argument 
implicates the expectation that the degeneracy or near degeneracy is lifted in second order of perturbation 
theory. However this is extremely implausible for thousands of nearly degenerate states. 

A generalization of scheme (ii) actually exists which extends its validity to all orders of the perturbation 
theory. This generalization, introduced by Brandow in a review article [12], is a Brillouin-Wigner (BW) 
type of perturbation theory, which requires to calculate the system energies self-consistently lfl3l . This 
self-consistency implies that BW expansions do not scale properly with the particle number in any finite 
order which can impede their use in solid state physics. 

A considerable number of perturbation techniques for nearly degenerate states were developed for pur- 
poses specific to quantum chemistry. Here we focus on a perturbation expansion which was developed 
by Lindgren lfT4l . It is of Rayleigh-Schrodinger (RS) type and can handle multi-particle problems with 
nearly degenerate ground states. From the formal structure it is more demanding than a BW expansion but 
it obeys the linked cluster theorem fl31 and the correct particle number scaling without adjustments. 

To our knowledge, neither the BW-scheme of Brandow nor the RS -scheme of Lindgren were used 
extensively in solid state physics in order to solve electronic problems with nearly degenerate states. In 
this context, it is of particular interest if these schemes can be employed effectively to investigate magnetic 
states of strongly correlated electron systems. Especially, can we reliably gain magnetic coupling constants 
if there is a competition between ferromagnetic and antiferromagnetic correlations and only higher orders 
in perturbation theory can settle the sign of the coupling? 

In the following section (Sec. [2}, we first introduce the concept of "nearly degenerate perturbation 
theory" and then present the schemes of Brandow and Lindgren for the generalization of BW and RS 
perturbation theory, respectively. As already mentioned, these schemes were little used in solid state 
physics and it is worthwhile to review the basics. In order to reassess their applicability in this context, 
we employ in Sec. [3] the different perturbation expansions to a two-site cluster with two orbitals and the 
standard set of local interactions, kinetic transfer of electrons, amended by a crystal field term acting in 
the two-orbital subspace. Such a term arises due to small deviations from the cubic symmetry, generically 
observed in the perovskites of transition metal oxides. In the Fock space with two electrons, this system 
may display a transition from a singlet to a triplet state for intermediate to large on-site Coulomb interaction 
U . This transition to a spin-orbital entangled state is not expected from conventional reasoning as it cannot 
be identified in standard second-order perturbation theory Ifl6l . We investigate the quality of the different 
perturbation schemes by comparing the 4 th order result to the exact solution. We summarize our analysis 
in the last section. 
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2 Formalism 

We consider the time-independent Hamiltonian problem 

H\* n ) = £U*„> 



(1) 



where H is in the Schrodinger representation. The eigenvalues E n and the eigenfunctions ^ n are labeled 
by n. We assume that the Hamiltonian H can be separated into 



H = H + H 1 , 



(2) 



where Hq is the unperturbed Hamiltonian which defines an exactly solvable problem and H\ is the pertur- 
bation. Accordingly, we assume that Hq can be diagonalized, and consider the energy corrections to the 
eigenenergy E^ corresponding to the eigenstate \n) of Hq. The eigenenergy with corrections up to an 

(k) 

appropriate order k will be denoted as En and, to infinite order, E n represents the exact eigenenergy of 
H. 
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Fig. 1 (Color online) Structure of the Hilbert space representative for the almost degenerate perturbation theory. It 
is comprised of a low-energy subspace D and a high-energy sector H/D. They are separated by an energy gap which 
is sufficiently large to allow a perturbative expansion in an "interaction" which relates states in D and H/D. The 
projection operators V and Q allow to project general state vectors onto D and H/D, respectively. The state indexed 
by n is the reference state for which we consider the perturbation expansion; it is always in D. 



Moreover, we suppose that the eigenstates of both, the total Hamiltonian H and the unperturbed Hamil- 
tonian Ho, can be assigned to either of the two sectors of the Hilbert space H: a low energy sector D of 
dimension d and a high energy sector H/D which is the orthogonal space of D (see Fig. Q]). We note, that 
the subspace D is often referred to as the "model space" (see, e.g., Ref. J7]). We assume that the states in D 
are all degenerate or nearly degenerate — with near degeneracy implying that the energy difference of the 
states in D is smaller than the difference to any of the states in H/B. This decomposition allows to devise 
a well-defined perturbation expansion if the perturbative part H\ of the Hamiltonian has matrix elements 
smaller than the gap between low energy and high energy states. 

The projection operators which are associated with the two complementary subspaces are: 

?> = E Ip><pI 

pen 

Q= E k)(?l = i-^ (3) 

geH/D 

where the index p (q) refers to orthonormalized eigenstates \p) (\q)) of Hq in the low (high) energy 
sector of the Hilbert space. The projections of the eigenvectors of H into the low energy sector are denoted 
by = 'P|\f r p). This notation will be kept for the rest of the paper. We use the normalization with 
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2. 1 Third order perturbation theory for nearly degenerate s tates 
We rewrite the Hamiltonian: 

H = H + (V + Q)Hi(P + Q)=H +VH 1 V + QH 1 V + VH 1 Q + QH 1 Q. (4) 

At this point we will simply redefine the unperturbed part of the Hamiltonian and the perturbation. With 
the definitions: 

Ho = Ho+VHiV, (5) 
Hi = VHiQ + QHiV + QHiQ, 

the perturbative problem is restated in an adequate fashion for nearly degenerate cases (6): Ho is to be 
diagonalized in the d dimensional subspace D and Hi constitutes the perturbation of these "zero order" 
states with energies £^ = (p\Ha\p). This reformulation of the problem has an irrefutably beneficial 
effect: all matrix elements of Hi connecting two states in D are zero by construction; therefore only states 
in the complement H/B can be intermediate states in the perturbation expansion up to third order. We 
apply the textbook formulae of non-degenerate perturbation theory for the new perturbation Hamiltonian 
"Hi to find for a state b£D 



njnqn/qn njnq lq _,q 1 q 2 njq2n 

pO) _ c(o) _ v th th i v 1 1 1 (fa 

t„ t„ -2. A£ + 2. A£ qi A£ q2 {b) 

q ' 91,92 Hi H2 

where we use the short hand notation A£ q = £n^ — £ q °^ and H^ q2 = {qi\Hi\q2), and 91,(72 are 
indices for states in H/D, exclusively. The first order energy correction (n\Hi \n) vanishes identically by 
construction because the state n belongs to the subspace D. The first order is already encoded in the setup 
of Ho, and the suppression of the first order in the subsequent perturbation expansion also simplifies the 
third order term to the form presented in Eq. (|6}. Beyond third order, intermediate states peD intrude and 
spoil the expansion due to the smallness of A£ p in the denominators. 

In order to overcome this difficulty, one has to recast the perturbation expansion into a form that allows 
to separate the virtual excitations into the high energy states from the low energy processes in such a way 
that virtual excitations into intermediate low energy states are excluded. Such a reformulation exists for 
the Brillouin-Wigner perturbation theory as well as for the Rayleigh-Schrodinger perturbation theory. For 
the first (BW) we refer to Brandow's approach iflZl and for the second (RS) to Lindgren's treatment lfT4l . 



2.2 Brandow's generalized Brillouin-Wigner perturbation theory 

The Brillouin-Wigner perturbation theory is straightforwardly generalized to account for a low energy 
subspace D of dimension d with nearly degenerate or exactly degenerate states: the states of the model 
space D are excluded from the intermediate summations in the perturbation series. In fact, this procedure 
is an exact resummation as shown in Sec. II of Brandow's review Ifl2ll and in Ref. f\M . For convenience 
we label this procedure by BBW (Brandow Brillouin-Wigner) as Brandow seems to have been the first to 
make this approach explicit for a "quasi-degenerate" model space (see also Refs. l9l [T7l[T8l ). 
In the standard BW approach the energy eigenvalues in D are found self-consistently from 

rrnq rjqn rr"9l rrtl92 rr92« 

E n = E^+Hr+ £ Hl \ + £ 1 (O) 1 1 (0) +••• CD 

i — ' TP _ TP^> — ' (TP _ Tp">\(TP — fi 1 ' >\ 
9 GH/D ^ n ^1 9i,9 2 eH/D V^n a qi JK-^n -^92 ) 

if D consists of a single state. Here, we use the short hand notation Hf iq ' 2 = (qi\Hi\q2). 

For the BBW scheme with a e?-dimensional model space D, an effective Hamiltonian lfT4l is introduced 
for each eigenvalue E n 

H£ s = VH T + VW n (8) 
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where the effective interaction V W n in the low energy sector is found from 



w n = h 1 v + h 1 r n w 



n 



(9) 



with the resolvent |[T9l 



\q)(q\ 



(10) 




which is exclusively defined in the high energy sector H/B. Eq. © allows to write the eigenvalue equations 



We emphasize that the states \q) are the unperturbed states of the complementary space H/B and they 
specify the effective Hamiltonian through the sum in T™ over the high-energy sector. The E n denote the 
perturbed eigenvalues in the low-energy sector which is spanned by the corresponding eigenvectors | \I>™) . 
The perturbed eigenstates j^Jj 1 ) are identified through 



Again, we stress that the resolvent T" and the effective Hamiltonian H™ s depends on the perturbed energies 
E n : There is an effective Hamiltonian H™ s for each E n , and the eigenvalue equation Eq. (fTTT > has to be 
solved self-consistently for each E n with the corresponding H™ s . 

As the states of the model space B are excluded from the intermediate summations in the perturbation 
series, it is obvious that infinitely many terms are missing in the BBW summations with respect to the stan- 
dard BW approach. However their respective contribution is generated through the initial diagonalization 
in the enlarged model space B. If B consists of a single state, BBW reduces to Brillouin-Wigner perturba- 
tion theory and if B covers the entire Hilbert space H, the diagonalization already solves the problem and 
the perturbation expansion collapses to the zeroth order term. In sec. [3] we will present a model where an 
appropriately chosen B reduces the full perturbation expansion to the second order term. In such a case 
no further higher order processes can be produced which relate states in the high energy sector H/B: due 
to the restriction to states in H/B on the summations, the second order term already induces the exact 
solution. 

Albeit the simplicity of the Brillouin-Wigner and BBW schemes, a formal disadvantage has to be faced: 
The linked cluster theorem does not apply [Q~2]|, which implies an incorrect scaling with the number of 
particles term by term in the perturbation theory. This inconsistency is only lifted in infinite order. 



Rayleigh-Schrodinger perturbation theory cannot be generalized for a problem with almost degenerate 
states in the model space B by excluding straightforwardly the states of B from the intermediate summa- 
tions in the perturbation series. Such an approach would be correct if the d states of B were completely 
degenerate. However we want to address the case with d almost degenerate states, where the degenerate 
subspaces of B may have the dimensions d a with 1 < d a < d and J2 a d a = d. 

In this situation it is convenient to introduce the wave operator fl [20 21 1, which transforms all unper- 
turbed states in B into the exact eigenfunctions of H, and expand £1 order by order 11141 . This yields the set 
of the desired d eigenenergies and eigenstates l22l to the appropriate order. As the approach traces back 
to the Rayleigh-Schrodinger perturbation theory, it does not suffer from a lack of extensivity, in contrast to 
Brillouin-Wigner perturbation theory. 

One introduces the wave operator il through 



(11) 




(1+ T n W n )\n) 

( i + r n H 1 + r n B x r n H x + r n H l T n H l T n H 1 + • • • ) \n) 



(12) 



2.3 Lindgren's nearly degenerate Rayleigh-Schrodingerperturbation theory 



\%} = n\*° n ) 



(13) 
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satisfying |* p ) = Q,V\^f p ) and 

nQ = 0. (14) 

Namely, if one projects any of the d low energy exact eigenstates | fy p ) of the full Hamiltonian onto the low 
energy Hilbert space D, the wave operator takes care of restoring the "out-projected" part of the eigenstate, 
while £1 yields a null result when operating on the complementary subspace. 
To zeroth order the wave operator is simply the projection operator onto D: 

n m = v (15) 

Once fl^ is obtained to order i — 1, the effective Hamiltonian H^g acting on D takes the form lfl4l : 

H { J = TH Q V + THi (n w + + fi (2) + . . . + fiC- 1 )) (16) 

The solution of the eigenvalue problem 

H^)=E^) (17) 

(i) 

yields the exact eigenenergies E n (to order i) in the low energy space, even though it operates on the 
eigenstates |n) of Hq in D. As the wave operator is unique for the entire low-energy Hilbert space, also 
the effective Hamiltonian is independent of the state n, in contrast to that of the Brillouin-Wigner type 
Hamiltonian in Sec. 12.21 

Starting from Schrodinger's equation, Lindgren obtained a recursion formula for fl^: 

i-i 

fl il) ,Ho\ = QH^ 1 - 1 ^ - fi ( '~ m) i?in (m_1) (18) 



The lowest orders are then explicitly obtained as: 



p(0) _ F (o) 



^ (3) = E Ei«><pi 



p(o) _ p (0) 
geH/DpeD 

am _ E Ew(p| Wflin ».- n c. Bin w- n m flin ».- nmgl |p) (19) 

<j£H/B pSB ^9 

With the knowledge of fiW, up to the appropriate order i — 1, one identifies the effective Hamiltonian if^£ , 
Eq. ( TToT ), and determines the eigenvalues En in the low-energy sector B through Eq. (fTTT i. 

Whereas the conventional approach of Sec. 12. H is a straightforward Rayleigh-Schrodinger perturbation 
expansion of the eigenenergies to the appropriate order in "Hi (see Eqs. © and ©), the Lindgren approach 
is an expansion of the effective Hamiltonian HQ to order i in H±. The diagonalization of H^2 then 

(i) 

introduces higher orders beyond i into the respective eigenenergies E n ' . 
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3 A two-site model to test perturbation theories 

It is now compelling to put the theoretical presentation into practical terms. A comparative analysis of the 
different perturbational schemes is realized for a basic model for which an exact solution is known: the 
two-band Hubbard Hamiltonian for interacting e g electrons on a two-site molecule. Such a model may be 
written as follows 



H — i?kin + -ffint + Hqy ■ (20) 

For e g electrons there are two orbital flavors: \x) ~ x 2 — y 2 and \z) ~ 3z 2 — r 2 forming a basis in the 
orbital space. Accordingly, the kinetic energy for hopping in y-direction is parameterized by 

# ki „ = £ ^'(c\ pa C 2p , a + h.C.), t""' =-{(^ f ) , (2D 

pp'a 

where t is an effective {(Ida) hopping matrix element, p is a band index with entries p = x or z, while a 
labels the spin states. The electron-electron interactions are described by the on-site terms, which we write 
in the following form 11231 

#int = U^(n tx ^n lxi + n^n lzi ) + (U - -J H )^2n lx n lz 

i i 

- 2J H ^ Six ■ &iz + 4)^ ( c ix1- C ixl C izl C iz1- + c iz1- C izl C ixj c ixl) ■ ( 22 ) 
i i 

Here U and Jji denote the intra-orbital Coulomb and Hund's exchange elements, whereas rii p = rii pa 
is the electron density at site i in the p = x,z orbital state. The term Hqf represents the uniform crystal- 
field splitting 

H C f = ^EcF^2{n iX a - n lza ) (23) 

resulting from, for example, a uni-axial pressure along the z-axis. 

Considerable interest in models with orbital degeneracy has been fueled by the investigation of the 
Mott transitions and their understanding (see, for example, Ref. l24l ). As concerns the lattice version of 
this particular model, the surprising richness of its phase diagram is instructive, especially in connection 
with the CMR manganites (for a review see l25l ). For instance, several antiferromagnetic phases and 
ferromagnetic phases have been studied in mean-field approximation on the square lattice ll26l . or by exact 
diagonalization Il27l and DMFT studies l28ll . The model is also believed to be a minimal model for the 
striped nickelates H29II30I . and striped phases have been studied in a number of ways (e. g. Refs. 01II32I ). 

In the following we consider 

Hi = H kin (24) 
as a perturbation of the local Hamiltonian 

Ho = H int + #cf- (25) 

With the proper perturbational expressions for the energy corrections, we can investigate the conditions 
under which the ground state of the model (f20) is ferromagnetic. In order to study this transition within 
the introduced perturbational schemes, one first needs to determine the set of states forming the subspace 
D. There are four subspaces, corresponding to states with either spin (one subspace), or spin 1 (three 
subspaces). Let us first focus on the triplet states, restricting ourselves to the ones with S z = 1. A basis 
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Fig. 2 (Color online) Energy of the lowest singlet and triplet states of the two-band Hubbard model for 2 sites with 
Jh/U = 0.2, Ecp/t = 0.5. The results for the Lindgren Rayleigh-Schrodinger scheme were generated with four 
nearly degenerate states in the low-energy subspace D. The upper two dashed curves display the result to fourth order, 
the lower two to second order. 



of the four-dimensional subspace D is given in Appendix A, and it is easily verified that these basis states 
are eigenstates of Hq. A basis of the two-dimensional complementary subspace H/D is also given in 
Appendix A. These states are eigenstates of Hq as well. We form the 2x4 matrix (q\Hi\p). It reads 



{l\Hi\v) 













V6 



(26) 



One can then easily evaluate Eq. (TT9b and solve Eq. ( fTTI ). If we now consider the singlet states, the calcu- 
lation runs in an analogous fashion. The subspace D is again four-dimensional, while the subspace H/D is 
now six-dimensional. Using the basis given in Appendix B, we obtain the 6x4 matrix (q\Hi\p) as: 



(27) 





( 


-V2 
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Fig. 3 (Color online) Results for the standard Rayleigh-Schrodinger and Brillouin-Wigner perturbative expansions 
for the two-band Hubbard model for 2 sites with Jh /U = 0.2, Ecf /t = 0.5. These evaluations use a one-dimensional 
low-energy subspace. 



Note that the states \q), which enter Eq. d27l i, are not eigenstates of Ho, and one further diagonal- 
ization step is necessary in order to obtain the desired (q\Hi\p) = 5Zm=i (qI^i l^ 1 ) i^ 1 \p)- This 
is described in Appendix B. Once these steps are worked out, the desired energy for eigenstate n to 
perturbative order i, viz. En \n), can be evaluated. The higher order terms of the wave operator in 
the Lindgren Rayleigh-Schrodinger technique simplify substantially: fl^ = = and = 

- EqGH/D E p gd \0)(p\ (<7|^ (1) -Hi^ (1) b)/(£p 0) - Ef ] ) due to the suppression of high-energy hopping 
processes. 

The result of this evaluation for both singlet and triplet states is shown in Fig. [2] Remarkably, the 
transition from the large U ferromagnetic to the intermediate U antiferromagnetic ground state that occurs 
at U c /t = 9.6 in the exact diagonalization calculation lfl6l (see vertical arrow in Fig. |2]i is reproduced 
in the perturbative treatment, i.e., in Lindgren's Rayleigh-Schrodinger type perturbation theory. It occurs 
here to fourth order perturbation theory at jyLmdgreiy^ _ g yg f QI f Qm near jy degenerate states in the low 
energy sector D. These states in D are identified by Eqs. d30b and (f28b for the singlet and triplet spaces, 
respectively. For second order perturbation theory with four nearly degenerate states [/Lmdgreiy^ = g o, 
(see Fig. |2]). The success of the perturbative evaluation derives from the correct identification of the (four- 
dimensional) low-energy space and the inclusion of the fourth order term, i.e. Q^ 3 \ 

Within non-degenerate Brillouin-Wigner and Rayleigh-Schrodinger perturbation theory there seems to 
be no transition for the chosen parameter set (see Fig. |3]for fourth order perturbation theory). As expected, 
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Fig. 4 (Color online) Energy of the lowest triplet state of the 2-band Hubbard model for 2 sites with Jh/U = 0.2, 
Ecf /t = 0.5. The Lindgren approach uses a 4-dimensional low-energy subspace D. In contrast, the displayed results 
for the standard Rayleigh-Schrodinger and Brillouin-Wigner perturbative expansions build on a one-dimensional D. 
The three perturbative approaches have been performed up to 4th order. 



these results are inferior due to the inappropriately chosen (one-dimensional) subspace D (cf. to the four- 
state Lindgren scheme and the exact result in Figs. [4] and |5). An exception is the evaluation of the triplet 
state within non-degenerate Brillouin-Wigner. 

The conventional approach J6), which was briefly presented in Sec. 12.11 reduces to the Rayleigh- 
Schrodinger perturbation theory for this specific model, as Ho = Ho in Eq. ((5]): there are no states in 
the low-energy sector which may be hybridized through VH\ V. Consequently, the conventional approach 
compares poorly with the Lindgren approach, even in second order (see. Fig. |6). 

The implementation of degenerate perturbation theory, for example, by setting Ecf to zero and perturb- 
ing the four-fold singlet state with respect to U and Ecf, which then lifts the degeneracy, has often been 
suggested. However, the artificial introduction of the zeroth-order degeneracy does not lead to adequate 
results, as seen in Fig.[5]for the evaluation up to fourth order. 

As the matrix elements obey (q\Hi\q') = for the considered model, the Brandow Brillouin-Wigner 
technique becomes exact for the four nearly degenerate states in the low energy sector D (given by Eqs. d30b 
and d28l )) already in second order. This is a special property of the chosen model and not generic for larger 
systems. 
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Fig. 5 (Color online) Energy of the lowest singlet state of the 2-band Hubbard model for 2 sites with Jh/U = 0.2, 
Ecf jt = 0.5. Both, the Lindgren approach and the degenerate perturbation theory, take a low-energy subspace D of 
4 states. In contrast, the displayed results for the standard Rayleigh-Schrodinger and Brillouin-Wigner perturbative 
expansions build on a one-dimensional D. All perturbative approaches have been performed up to 4th order. 



4 Summary 

Perturbation theory for atomic clusters is, for typically thousands of low energy states, genetically a per- 
turbation expansion with a large number of almost degenerate states. Moreover, third or fourth order 
evaluations are often necessary to attain the required accuracy or include new qualitative effects such as 
electronic ring exchanges. Consequently, a perturbation-theory scheme for almost degenerate states is an 
essential tool in theoretical physics. Almost degenerate perturbation theory has been elaborately investi- 
gated in quantum chemistry but, to our knowledge, these approaches have had little impact on solid state 
modeling. In fact, second or third order theory is well established (6), but it cannot be straightforwardly 
extended to fourth order and beyond (see Sec. 12. I) due to intruding low-energy states as intermediate 
states. A much cited technique is to start with the corresponding degenerate problem and then introduce 
the deviation from the degenerate case as perturbation; such an approach could benefit from the commonly 
accepted schemes for highly degenerate perturbation theory. However, this approach leads to unsatisfying 
results in higher order, on account of systematically inaccurate energy denominators — this shortcoming is 
confirmed for the investigated model (see the result for "degenerate perturbation theory" in Fig. |5}. 
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Fig. 6 (Color online) Energy of the lowest singlet state of the 2-band Hubbard model for 2 sites with Jh/U = 0.2, 
Ecf /t = 0.5. Both, the Lindgren approach and the conventional perturbative expansion J6] have been performed up 
to 2nd order. 



A suitable perturbation theory has to start from unperturbed low- and high-energy states which are 
diagonalized in the respective subspaces. One then has to devise a scheme in which the low-energy sub- 
space is not retraced in the "string" of virtual excitation processes induced by the perturbation, otherwise 
small energy denominators spoil the evaluation. This is achieved in Lindgren's scheme lfT4l for a per- 
turbation theory of Rayleigh-Schrodinger type (see Sec. 12.31 ). In this scheme, an effective Hamiltonian 
is to be identified through appropriate traces over high-energy states, and the eigenvalue problem for the 
effective Hamiltonian provides the perturbed eigenenergies and states in the low-energy space. For a 
Brillouin-Wigner type perturbation theory the approach is methodically very simple: the low energy in- 
termediate states are excluded in the perturbation-theory summations, as originally presented by Brandow 
(see Sec. 12.21 ). The Lindgren scheme can be derived from the Brandow Brillouin-Wigner expressions 
through expansion j 12). The schemes were introduced in nuclear physics and quantum chemistry but have 
not received much attention in solid state physics. 

For the basic example of a two-site two-orbital Hubbard model we presented a comparative analysis 
of the most prominent perturbative approaches and examined them against a singlet-triplet transition for 
strong electronic correlations which is known from the exact diagonalization. The Lindgren approach 
reproduces this transition even quantitatively (in fourth order). The Brandow scheme becomes exact in 
second order if the subspace of low-energy states D is appropriately chosen. However, this property of 
the Brandow perturbation expansion is particular for the considered model and will generally not persist 
for larger systems. As to the formal structure, the Lindgren Rayleigh-Schrodinger expansion is more 
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demanding but it obeys the linked cluster theorem in each order of the expansion which implies the correct 
particle-number scaling. 

Acknowledgments. R.F. gratefully acknowledges the Region Basse-Normandie and the Ministere de 
la Recherche for financial support. This work was supported by the Deutsche Forschungsgemeinschaft 
through TRR 80 (T.K.). C.H. and T.K. are grateful for the warm hospitality at the CRISMAT in Caen 
where part of this work has been done. 



A Basis of D in the triplet subspace 

When considering triplet states (with S z = 1), a convenient basis for the subspace D is given by {\p) T , 1 < 
p < 4}, with : 



|1)t = 


c i,z,T c 2,2,T 0) 


|2) T = 


— ( r f r f 


|3) T = 


— ( r f r f 


|4> T = 





C l,z,t C 2,a;,t) I") 
C l,2,t C 2,a:,t) 1^) 



(28) 



These states are clearly eigenstates of Hq, with eigenvalues i?j ' 1 = —Ecf, = = 0, and 

E^ = Eqf- As for the complementary space H/D, we use the basis {\q) T , 1 < q < 2} given by: 

i /j )t = ^(4,^,t-4,,,t c L,t)i >- ( 29 ) 

Both are eigenstates of H , with eigenvalues U — 2>Jh- 



B Basis of D in the singlet subspace 

When considering singlet states, a convenient basis for the subspace D is given by {\p) s , 1 < p < 4}, 
with : 



!>S = 4( c l z ,T c L 4 - ct M4 c L,t)l ) 



V2 

|2)g = ( r 'l ..< .-' 2.-... — C l,a;,| C 2,z,t C l,z,T C 2^4 _ C \ ,z ,i°2 ,x ,\) 1^) 



l^)s — 9 ( C l,x,t C 2,z4 C l,a;,| C 2,z,t C L,t C 2,i,4 + C l,z^ C 2,a;,t) 1^) 



l 4 >5 = ^( c U t 4,^-4,x4 c Lt)|0>- (30) 
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These states are eigenstates of Hq, with eigenvalues E\ ' = —Ecf, E. 2 ' = E$ = 0, and E\ ' = 
Ecf- As for the complementary space H/D, we start with the basis {\q) s , 1 < q < 6} given by: 



Unfortunately these eigenstates of H- m t are not eigenstates of Hq, as 5 (I\Hq \II) s = s(IH\Hq\IV) s = 
—Ecf, but the resulting 2x2 blocs are most easily diagonalized. As a result we obtain the eigenvalues of 
the unperturbed Hamiltonian as e\ 0) = e\% = U - ^ J% + E' CF , E*f} = E<$ = U + ^ J' 2 H + E CF , 
andE^ = = U - J H . 
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